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Abstract 

We give a rigorous treatment to the thermodynamic limit of the 
6-vertex model. We prove that the unique solution of the Bethe- 
Ansatz equation exists and the distribution of the roots converges to 
a continuum measure. We solve this problem for < A < 1 using 
convexity arguments and for large negative A using the Fixed Point 
Theory of appropriately defined contracting operator. 



1 The 6-vertex model and formulation of the 
problem 

The 6-vertex model is an exactly soluble model of classical statistical me- 
chanics introduced and solved in various special cases by Lieb [U EJ A 
solution of the most general case was obtained by Sutherland |4J. A clear 
description of this model and various other soluble models can be found in 
Baxter's book [3]. However, as Baxter remarks, an exact solution is not the 
same as a rigorous solution. In fact, already in his first article on the ice 
model pp, Lieb initiated the rigorous analysis of the model. A more exten- 
sive analysis was made by Lieb and Wu [6J. An important technical question 
was left unresolved, however. This concerns the convergence of the distribu- 
tion of (quasi-) wavenumbers to a continuum measure in the thermodynamic 
limit. (Another technical problem, i.e. the independence of the free energy 
on the boundary conditions, was resolved by Brascamp et al. [7j.) A similar 
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problem was solved in the case of the nonlinear Schroedinger model in [8]. 
The 6-vertex model is more complicated because we cannot in all cases use 
the convexity argument of Yang and Yang [0] used there. However, their 
argument does extend to a certain domain of parameter space. Here we 
show how it can be used to prove the convergence of the Bethe Ansatz solu- 
tions in the thermodynamic limit in that case. In addition, we use another 
technique for proving the existence of a unique solution to the Bethe Ansatz 
equations in the thermodynamic limit in a different domain of parameter 
space. Uniqueness in other parts of parameter space is still an open problem, 
though numerical iteration does seem to converge to a unique solution. 



1.1 Definition of the model and the free energy 



We first recall the definition of the 6-vertex model and some general re- 
sults concerning the existence of the thermodynamic limit. We then review 
the transfer matrix formulation of the model and the diagonalisation of the 
transfer matrix by means of the Bethe Ansatz. 

The 6-vertex model is a model of classical statistical mechanics where the 
configurations are given by arrows on the bonds of a 2-dimensional square 
lattice. At each vertex only six different configurations of arrows are allowed 
(the so-called ice condition): 



e l e 2 e 3 

Each of these vertex configurations is assigned an energy and we assume 
spin-flip invariance, so that the first and the second, the third and the fourth 
and the fifth and the sixth configuration have the same energy. We denote 
these energies by e%, e 2 and e 3 . If (3 is the inverse temperature, the cor- 
responding Boltzmann weights are: a = exp[— /3ei\, b = exp[— /3e 2 ] and 
c = exp[— (3e%\. The partition function is therefore 



Z MtN (a,b, c) 



reCjvi.iv 



where M is the number of rows and iV is the number of columns in the 
lattice, Cm,n denotes the set of allowed configurations, and the total energy 
of a configuration T £ Cm,n is 

EiY) = m(r) e 1 + n 2 {V) e 2 + n 3 (V) e 3 , (1.2) 
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if n,(r) is the number of vertices of type i in the configuration 17 "Solving 
this model" now means: finding an explicit expression for the thermodynamic 
limit of the free energy density, i.e. 

f(e 1 , e 2 , e 3 ; (3) = -- ^hm^ — In Z N , M {a, b, c) (1.3) 

We shall assume periodic boundary conditions. 

The first question that arises is whether the limit (11.31) exists. This was 
solved by Lieb and Wu [6J. In fact, we need the existence of the free energy 
at constant density p. In [6] this is highlighted as an open problem, but in 
fact, in the case of periodic boundary conditions, their method extends to 
this case. For convenience we repeat their argument here. (It was proved 
by Brascamp et.al. [7j that periodic boundary conditions are equivalent to 
free boundary conditions in the thermodynamic limit. It should be noted 
that not all boundary conditions are equivalent: see the recent solution of 
the model with domain-wall boundary conditions by Bleher et al. [TO].) 

The periodic boundary conditions imply that, in a given configuration, 
the number of up arrows in every row of vertical arrows is the same. We 
shall call this number divided by the maximum number N, the density p. 
The partition function with fixed density p is given by 

z M , N (w 1 ,w 2 ,w 3]P )= E E E e-mr) ' ( L4 ) 

a:#{ ft ,=+i}=jvp 7 ree p M>N (a,j) 

where C M N (a,j) = Cm,n(&, 7, 7) and Cm^n(&, ot' , 7, 7') denotes the set of 
configurations with given boundary arrows: a and a' for the bottom and top 
rows of vertical arrows, and 7 and 7' for the left- and right-hand columns of 
horizontal arrows. 

Proposition 1.1 Let Z P M N denote the partition function of the six-vertex 
model with periodic boundary conditions and let (Mi,Ni) be a sequence tend- 
ing to infinity in the sense of Van Hove, and suppose that (pz)£i ^ s a sequence 
of numbers pi G [0, 1] tending to p such that piNi e N. Then the correspond- 
ing free energy density f p (ei, e 2 , e 3 ; (3, p) defined by 

f p {e 1 ,e 2 ,e 3 ;f3,p) = -- lim — ^— In Z P M N (a, b, c; pi) (1.5) 

p l^oo MiTsi 

exists and is independent of the sequences (Mj, iVj) and (pi). Moreover, 
f p (ti, £2, £3; P, p) is convex as a function of p and concave as a function of 
the variables e±, e 2 , 63 and (5. 
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Proof. We start by considering special sequences. Assume first that p G 
[0, 1] n Q. Take N G N so large that pN G N, and choose M G N arbitrary. 
Consider a standard sequence of rectangular boxes of height Mi = 2 l M and 
width Ni = 2 1 Nq. One then proves as in Lieb and Wu [6] that the limit (II. 
exists, using the inequalities 

Zm 1> n 1 (P,p)>&m 1 j* 1 {.P,p) > (Ojui-Lty-iCAp)) 

> (^M i _ 1 ,iv i _ 1 (/5,p))2- 4 ^ + ^). (1.6) 

(We suppress the dependence on ej and on the periodic boundary conditions.) 
Here we define 

Qmn(I3,p)= max max V e _/3JS(r) . (1.7) 
The inequalities (11 .6I> imply that the sequence 

fMuNA^P) = ~ 0M lNl lnZ Mi,N,(P>P) ( L8 ) 

is essentially decreasing. As it is also bounded below, it converges. The 
concavity as a function of e$ is standard. To prove the convexity as a function 
of the density, suppose that p\ < pi- Then if ol\ and a 2 are given sets of N 
vertical arrows with densities p\ and p% respectively, and we write a = «iUo2 
for the union, 

£ Y, e~ mT) >J2 S Yl e -^ E ^ +E ^\ (1.9) 

T reC M,2iv( Q .7) 7 r 1 ecP fiV (ai,7)r 2 eC^ iJV (o i2l 7) 

Summing over a x and we have 

+ p 2 )) > ^,Ar(Pi)^M,Ar(P2) (1-10) 

and hence 

/(/3, |( Pl + p 2 )) < i (/09, Pl ) + /(/?, P2 )) . (1.11) 

Convexity implies continuity and we can thus extend the definition to all 
P G(0,1). 

To show that the definition of f(/3, p) is independent of Mo and iVo we fill 
a general domain A with rectangles and use the condition e\ > and e 2 > 
to decorate the remainder by vertices 1,2,3, and 4, as in [5]. □ 

In the exact solution of the six-vertex model one actually takes the limits 
M — > oo and N —>■ oo consecutively, but it was also shown by Lieb and Wu 
[6] that, for periodic boundary conditions, this yields the same limit as (11.51) : 
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Proposition 1.2 The double limit 

/(ei,e 2 ,e 3 ;/3,p) = -- lim lim ——hxZ^ N (a,b,c) (1.12) 

p N-^oo M^oo 1\ lvl ' 

exists and equals f(e 1 ,e 2 ,e 3 ,P,p). 

1.2 The transfer matrix and its diagonalisation 

The transfer matrix method for solving models of classical statistical mechan- 
ics is common knowledge. Using periodic boundary conditions one writes the 
partition function as a trace 

Z M , N = Trace (V^) (1.13) 

where Vn is the transfer matrix with entries between two rows of vertical 
arrows a and a' given by 

N 

(^W = EII ex p[-^(^'^)]- (i-i4) 

The sum runs over a row of horizontal arrows 7 = (71, ... , 7^), where 7„ is 
the horizontal arrow between the (n — l)-th and the n-th vertex. It follows 
that we can take the limit M — > 00 to obtain 

/CM = "4 j im ^lnA max (iV), (1.15) 

p AT— >oo iv 

where A max (iV) is the maximum eigenvalue of the transfer matrix Va? which 
exists because Vn satisfies the conditions of the Perron-Frobenius Theorem. 
The transfer matrix can be diagonalised by means of the Bethe Ansatz. If we 
write \xi, . . . , x n ) for the row configuration with n up-arrows then a general 
wave function in the subspace with n up-arrows can be expressed as 

ip = ^2 il>(xi,...,x n )\x 1 ,...,x n ). (1.16) 

l<xi<-<x„<JV 

The Bethe Ansatz for eigenfunctions of Vn then reads: 

n 

ij)(xi,...,x n ) = ^ A a exp [j fc^gj]. (1.17) 

Here, the sum runs over the set S n of all permutations of{l, . . . , n) and the 
coefficients A a and the wave numbers ki,...,k n are to be determined by 
inserting into the eigenvalue equation. This yields the following conditions: 
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1. The wave numbers must satisfy the simultaneous nonlinear equations: 



e iNkj 



-l)™- 1 [ [ e - ie ^ M \ (1.1? 



where the function 9 is defined by 

1 _ OA P ik _|_ 

exp [—i9(k, k')} = 1 _ 2Ae , tfc ,; e , fc+fc , ) (1.19) 

with 

a 2 + b 2 - c 2 . 

A = , . 1.20 

2ab v ; 

2. The corresponding eigenvalue is given by 

n n 

A(h, ...,k n ) = a N H L(e^) + b N J] M(e^), (1-21) 

3=1 3=1 

where L(z) and M(z) are given by 

a 2 — abz 
. « a 2 — c 2 — ate 

M < 2 > = a6 -^ • < L23 > 

Of course, (11.191) only defines the function 9 up to a multiple of 2tt. In 
taking the logarithm of (ll.lSp . we shall assume that — 7r < 9(k,k') < 7r. We 
obtain 

n 

iVfy = 2irlj - e ( k j, h), (1-24) 
i=i 

where Ij G Z if n is odd, and Ij G Z + | if n is even. These equations 
are identical to the BA equations found by Bethe [UJ in his solution of the 
Heisenberg chain. They were analysed in detail by Yang and Yang [9], who 
showed that the ground state of the Heisenberg chain is obtained by choosing 

I j= j-l(n + l). (1.25) 

They also showed that, for this choice, the equations (11.241) have a real so- 
lution for ki, . . . , k n . Lieb [1] then argued that, as the Heisenberg Hamil- 
tonian also satisfies the conditions for the Perron-Frobenius Theorem, the 
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corresponding eigenfunction must be positive, and hence it must also be the 
eigenfunction of the transfer matrix with maximum eigenvalue. We therefore 
have 

Amax = A(A;i,...,Av,) (1-26) 

where k%, . . . , k n are the solutions of (11.241) in case the Ij are given by (11.251) . 
By formula (11.151) . the free energy is now given by 

MP) = Jim min | Cl - _5>L( c *<),e a - _5>M(e<*) j ■ 

(1.27) 

By the fact that the free energy is convex and symmetric in the density, the 
minimum is attained at p — \ and the solution for k\, . . . , k n corresponding 
to the integers ( 11.2511 . 

In this paper we address the question of how to compute the thermody- 
namic limit (11.271) . We want to take the limit iV — > oo, keeping p = n/N 
fixed. One usually makes the reasonable assumption that, in this limit, the 
distribution of the wavenumbers ki, . . . , k n tends to a continuous distribution 
with density p(k). In the following we shall investigate the validity of this 
assumption. Following Yang and Yang |9], we consider separately the cases 
A e [0, 1) and A < 0. (The case A > 1 is trivial.) In the attractive case, 
A G [0, 1), we can apply the same reasoning as in the case of the nonlinear 
Schroedinger model (see [H]) and use the convexity of a certain functional to 
prove the existence of a unique solution to (11.241) . In the repulsive case, we 
can only treat the case A <C — 1. The case of smaller negative A is more 
delicate. The difficulty is proving the uniqueness of the solution. We can 
show, however, that if one assumes that the solution is monotone, the lim- 
iting solution is unique. Numerical solution of the BA equations seems to 
suggest that it is unique even without this assumption, but we have so far 
been unable to prove that. 

2 Thermodynamic limit in the case A £ [0, 1). 

In taking the thermodynamic limit we distinguish the cases A > 1, A e [0, 1), 
A e (-1,0) and A < -1. The case A > 1 is trivial (Cf. Baxter [5]) 
so we start with the case A G [0,1). We first prove an analogue of the 
existence and uniqueness of a solution to the Bethe Ansatz equations in the 
thermodynamic limit. In the present case this is analogous to the nonlinear 
Schrodinger problem treated in [8]. 
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Theorem 2.1 Let m E Ai b + [— §,§] with \\m\\ < 1/2 and supp(m) C 
[— ti\ \m\ |, 7r| \m\ |]. In case ||m|| = | ; assume that there exists 5q > stzc/i 
i/iai /or < 5 < 5 Q , 

-({ 9£ [-|,|]:|-M<*})<^ (2.28) 

(Notice that the uniform distribution satisfies this condition.) Let A = 
— cos with /j, e (7r/2,7r). T/ien i/iere exists a unique continuous function 
k : [— 7r/2, 7r/2] — > [— 7r + /i, 7r — /i] swc/i that 



r w/2 

k(q) = q- 6(k(q),k(q'))m(dq'). (2.29) 

J-n/2 



-tt/2 
t/ 

Proof. Define the new function g(g) by 

P */2 _ pfl(g) 

e *W = e f . (2.30) 

Then = K(g(q)) where K : R — > (— 7r/2,7r/2) is an increasing func- 
tion given by 

*(«) = f ^ = 2 tan- ( ^^13.) . (2 .3l) 

v ; J cosh(/3) -cos// V tan(///2) / v ; 

It follows that g(q) must satisfy: 

K(g(q)) = q ~ [ ' uj(g(q) - g{q')) m(dq') (2.32) 

J-w/2 

where 

, / tanh(a/2) \ 

-w=- 2t - (^r)- < 2 - 33 > 

Notice that 

,. , sin(2u) 

u/ (a) = : K - J - J - — r > 0. 2.34 

cosh a — cos(2/x) 

As in jH], we now define a functional B[g] on the space L 2 (IR, m) by 
B[g) = J S{g(q))m{dq)- J qg{q)m{dq) 

+ l -J j 'n(g(q)-g(q'))m(dq)m(dq'), (2.35) 
where S(a) = j Q a K{f3)d(3 and tt(a) = j* u{/3)d/3. 
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The functional B is well-defined because < S(a) < ^K'(0)a 2 and 
< fi(a) < \uj'{Q)a 2 } where 

K'(0) = , and o/(0) - 



1 — cos(/i) 1 — cos(2yu) 

It is also easily seen to be continuous. The Gateaux derivative in the direction 
of a function / is given by 

DB[g}f = j {#(</(?)) -g + J u(g(q) - g(q'))m(dq')j f(q)m(dq). (2.36) 

It follows that the solution to (3.4) is a stationary point of B. Moreover, B 
is convex as 



dt 2 



B[g + tf] = j K\a(q))f(q) 2 m(d q ) 



+1 I I u'(g(q) - g(q'))(f(q) ~ fW)) 2 ™(dq) m{dq') > 



2 




(2.37) 



by ( 12.34p and the fact that K'(a) > 0. This proves the uniqueness of the 
solution. To prove the existence, we need to find a compact set which contains 
the minimiser. 

Consider first the case that ||m|| < |. Now, as a — ► ±oo, K(a) — > 
±(7r — fi) and uj(ol) — > ±(2/i — 7r). Let M be so large that 7r — /i — |if (a)| < e 
and (2/i — 7r) — |o>(ck)| < e for |a| > M, where e > is to be determined later. 
Consider the set 



M 



{q E [-7r||m||,7r||m||] : g(q) > M}. (2.38) 



For M large enough, we can assume that m(YM) < e - We now replace g on 
the set T 2 m by ±2M, i.e. we set 

g(q) = sgn(g(q))mm{\g(q)\,2M}. (2.39) 
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By convexity of the functions £1 and S we then have 
B[g]-B[g] = 

= J (S(g(q)) ~ S{g(q))) m{dq) - J q {g{q) - g(q)) m{dq) 

+ l -J J (n(g(q) - g{q')) - Q(~g(q) - g(q'))) m(dq)m(dq') 

> I (\g(q)\-2M)(K(2M)-\q\)m(dq) 
J r 2 M 

+ I m(dq) j m(dq') (Q(g(q) - g(q')) - Q(g(q) - g(q'))) 

(2.40) 

where we used the convexity of the function S and the fact that if q, q' G Y 2 m 
then the second term in the double integral is zero whereas the first term is 
positive since Q > 0. Next using the convexity of fl and the above bounds 
on the derivatives we get 

B[ 9 ] - B[g] = 

> j (\g(q)\-2M)(K(2M)-\q\)m(dq) 
J r 2 M 

+ j m(dq) j m(dq')u(2M - | <?(?') I) (I <?(?)! - 2M) 

> j (\g(q)\-2M)(n-pi-\q\-6)m(dq) 
m(dq) [ m(dq') ((2/i - tt) - e) - 2M) 

> [ rn(dq)(\g(q)\ - 2M) (tt - pt - \q\ - e + ((2/i - tt) - e) m{V c M )) 
J r 2 M 

> j m(dq)(\g(q)\ - 2M) (tt - /i - \q\ - e + ((2/i - tt) - e) (||m|| - e)) 
" r 2 M 

> / rn(dq)(\g(q)\ - 2M) ((tt - /i)(l - 2||m||) - e(l + 2/i - tt + ||m||)) 
J r 2 M 

> (2.41) 
provided 

(7r-/i)(l-2||m||) 



e < 



1 + 2fi — 7i + | \m\\ 

We conclude that the minimiser must satisfy \\g\\oo < 2M and is a fortiori 
contained in the ball {g G L 2 (m) : ||g|| 2 < M}. This ball is bounded and 
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therefore weakly compact. But the functional B[g] is norm continuous and 
convex and therefore lower semicontinuous for the weak topology, see e.g. 
|12j . Prop. 1.5 of Chap. 2. It follows that it attains its minimum on a 
compact set. 

Next consider the case ||m|| = |. In that case we cannot prove that the 
minimiser is bounded, so we need a more sophisticated bound. We use the 
function 

/(g) = -2 In (|-|?|). 
Given M > and 5 > 0, we define the sets 



= {?€ hi*, I*] : > M, \q\ < \ - 6} (2.42) 

and 

T k = {qe [-Itt, §tt] : |s(g)| > /(g), | - ^ k+1 6 < \q\ < | - 7 ~\f} , 

(2.43) 

where 7 > 1 is a parameter to be determined later. 
We now write 

r M = if u |J r k 

k>l 

and consider the decomposition 

{(q,q') GhfTT, Itt] 2 : g G T M or g' G T M } = 

= U( rfcX (U r M u (U r <) xr fc u(r fc xr fe )J. (2.44) 

fc>0 V \l>k / \l>k / / 

Note that this is a disjoint union. Replacing now g{q) by 

g(q) = sgn(#(g)) min{|#(g)|, (f{q)xvM\ V M +2Mx f m)} 
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we have first of all 

(S(g(q)) - S(g(q))) m(dq) - J q (g(q) - g(q)) m(dq) 
> [ m{dq){\g{q)\-2M)(K(2M)-\q\) 



2M 




+ J2 m(dq)(\g(q)\-f(q))(K(f(q))-\q\) 

k=l jTk 

> [ m(dq)(\g(q)\-2M)(ir-n-\q\-r]) 



CO „ 

+ J2 / m(dq)(\g(q)\-f(q)) 
k=i ^ r fc 



^I'-'-ipitri' 1 )'- 1 '!!' (2 - 45) 



where we used the bound 

which follows from the inequalities 

tan _1 (x — 5) > tan _1 (x) — 5 

and 

tanh(:r) > 1 - 2e~ N . 

For the term 

\f J m 9 (q) - g{q')) - n(g(q) - g(q'))) m(dq)m(dq') 
we consider the contributions from the decomposition (I2.44p separately: 

(Q(g(q) - g(q')) - Q(~g(q) - g{q'))) m{dq)m{dq') 
> I m(dq)(2n-n - rj)m((T M ) c ) (2.47) 

as before. Combining this with the first term of (12.451) gives a positive con- 
tribution provided M is so large that m{T M ) < e and K{2M) > it — fi — r] 
and uj(M) > 2/i — n — r\ where \r\ + ne < 5. 
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Next consider a term of the form 

(n(g(q) - g(q')) - Q(g(q) - g(q'))) m(dq)m(dq'). 




'i* ^(u ; > fc r fe ) c 
Assuming 5 < 8 , this is bounded by 




r*; ^(U(>fer fc ) 
> 



x 



{M{g{q) - g{q')) - Sl{g{q) - g{q'))) m{dq)m{dq') 

m(dq) (\g(q)\ - f(q)) 
4 



2/i — 7T 



Since 



l>k 



tan(7r — fi) 
iC<qE 



??? 



\l>k 



(2.48) 



7T 7T 



2 ' 2 J 2 



7T 



7 - fc+1 5 < |g| 



we have by the assumption about m, 



rn 



(\Jr,)< V+v 



(2.49) 



Therefore 




(to(9(q) - 9(4')) - ^(g(q) - g(q'))) m(dq)m(dq') 



> 



x 



m{dq) {\g{q)\ - f(q)) 
4 



2n-«- j— r (--\ q \) 2 ](-- ~l- k+1 s) 

p tan( 7 r- / u) V2 m ) \ \2 tt ' J 



j- k+1 5 . (2.50) 



Combining this with the corresponding term of (12.451) we have 
(S(g{q)) ~ S(g(q))) m{dq) - J q {g{q) - g{q)) m{dq) 
+ I I (n(g(q) - g{q')) - {l(~g(q) - ~g{q'))) m(dq)m(dq') 




T fc ./(Ui> fc r fe ) c 
> / m(dq)(\g(q)\-f(q)) 



71 



X 7T — fl — 
+ [2(1-71 



7T 



tan(///2) V2 

4 (TX 



l?l - M 



tan(7r — /i) V 2 
13 



-Ifll 



1 1 

2 ~ ^ 



.(2.51) 



Since f — 7 < \q\ < f — 7 fc 5 for g G T fc , we have 



where 



Hence 



c - 4 c - 4 andc-c X c 

Cl tan(///2)' ° 2 tan(7r — //)' &n ° Cl 2° 2 ' 



(2.52) 



(S(g(q)) - S(g(q))) m(dq) - j q (g(q) - g(q)) m(dq) 

I I (tt(g(q) - 9 {q')) - n(~g(q) - g(q'))) m(dq)m(dq') 

> [ rn(dq)(\g(q)\-f(q)) 



X 



(2.53) 



Finally consider the terms 



^ m(dq) / m(dq')(n(g(q)-g(q'))-n(f(q)-f(q'))). 



Since < Q(a) < (2/i — 7r)|a|, these can be bounded by 



m(dq) / m(dq , )^(g(q)-g(q / ))-^(f(q)-f(q f ))) 
r fc Jv k 

>-(/m-1) f m(dq) f m(dq')\f(q)-f(q')\ 

> it ln( 7 - fe (5)m(r fe ) 2 > 7r( 7 - fe 5) 2 1x1^5). (2.54) 
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In all, we get 

(S(g(q)) - S(g(q)))m(dq) - I q (g(q) - g(q)) m(dq) 



+ / / (nfo(g) - g(q')) - n(g(q) - ~g{q'))) m(dq)m(dq') 

+\( m(dq) f m(dq')(Q(g(q)-g(q'))-n(f(q)-f(q'))) 



2 

> I m(dq)(\g(q)\-f(q)) 



k 



7T 



(2.55) 



Choosing 7 < 2fl ~ 7T (which is possible as // < tc) and 5 small enough, this is 
positive. 

It now follows that in this case the minimiser of B[g] must satisfy 

\g{q)\ < f(q)x[-i 7r ,-L 7r +8]u[^-8,^] + 2M X[-i 7r+5 ,in-6] 

and therefore 



\\g\H < 2M 2 + 2 f(qym(dq) 
8 f s 

< 2M 2 + - (lnx) 2 dx < +00. (2.56) 
K Jo 

Again, it follows that B[g\ attains its minimum on this compact set. 

We finally prove that the unique solution g G L 2 (R, m) of (I2.32p in fact 
has a continuous version as a function g : [— ir\ \m\ |, ir\ \m\ \] — > [—00,00]. We 
have shown that there exists g G C 2 (R,m) satisfying (12. 32ft for m-a.e. q. We 
define the image measure fh = g(m) and put 



h(x) = K(x) + j oo(x — x')m(dx). (2-57) 
Clearly, h is a C°°-function on M and 

h'(x) = K'(x) + / uj'(x - x')fn(dx) > 0. 



Therefore, the inverse function g = h^ 1 is well-defined and C°° on the range 
of h. Since 

h(x) — > ±(ir — ji + (2/i — 7r)||m||) 
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— > ±00, the function g is defined on the interval 

An = (-0 - A* + ( 2 M ~ 7r )ll m ll)> 7r - A* + - 7r)||m||). 
Notice that if ||m|| < |, this interval contains [— ir\ \m\ |, 7r|m| |], whereas if 



m 



|, J m = (— |7r, |7r). In the latter case, g extends continuously as 
a function g : [— §7r, ±7r] — >■ [—00, +00]. Inserting x = g(g) we have for 
q G supp (m), 



h{g{q)) = K(g(q)) + / co{g{q) - g{q')) m(dq') = q 
for m-a.e. q. Hence g(q) = g(q) for m-a.e. q. Now inserting x = g(q) we get 
q = h(g(q)) = K{g{q)) + J ' uj{g{q) - g{q'))m{dq') 
= K{g{q))+ uj{g{q)-g{q'))m{dq r ), 



so that g satisfies ( 12.321) for all q in its domain. It remains to show that 
the solution g is unique. It follows from the mean-value theorem that any 
continuous solution is differentiable and its derivative is given by 

9\q) 



K>{g{q)) + fu'{g{q)-g(q'))m{dq>) 
1 

K'{g{q)) + j u'{g{q) - ~g{q'))m{dq'Y (2 ' 58) 

Note that the function g is uniquely defined modulo an m-null-set, so that 
the right-hand side only depends on the value of g at q. Since g(q) is uniquely 
defined on supp (m) by continuity, its extension to [— n\ \m\ |, tt| |m| |] is also 
unique. 

□ 

Theorem 2.2 The mapping m 1— > k m defined by Ii2.29\) in Theorem \2.1\ is 
continuous, that is, if m n — > m weakly then k mn k m in norm. 

Proof. Let rrin be a subsequence. Notice that ||fc m „|| < tt — fi and k mn 
is also equicontinuous because 

± 9(k k >) = A cos ( fc/ ) + CQS M > 

dk v ' ; A 2 sin 2 (A;- k')/2+ [cos(fc + fc')/2 - Acos(A;- k>)/2] 2 ~ 

(2.59) 
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for —Ti + fj, < k' < Ti — ji. Hence 



k'mM = I 1 + / jfcV^b) - e (°- 1) ( 2 - 6 °) 

Therefore, \k mn (q) — k mn (q')\ < \q — q'\ uniformly in n. It follows that there 
exists a subsequence m$ of m$ such that k (t, converges to a continuous 
function k uniformly on [— 7r/2, 7r/2]. We must show that k = k m . But 
9 is uniformly continuous on [— n + /i, tt — /i] 2 so 9(k m (2)(q) — & m (2)(-)) — > 
9{k(q) — &(•)) in norm, and hence 



It follows that /c(g) = k m (q). 



9{k{q) - k{q'))m{dq'). 



□ 



Corollary 2.1 // m„ — > m weakly, and m n is the image measure of m n 
under the mapping k mn then rh n — > m = k rn {m). 

Proof. Let F e C([-ir + tt - /j]). Then / F(k)m n (dk) = 
I F(k mn (q))m n (dq) and 



J F(k)m n (dk) - J F(k)m(dk) < 
< J \F(k mn (q)) - F(k m (q))\m n (dq) 

F(k m {q))m n (dq) - / F (k m (q))m(dq) 



+ 



(2.61) 



The right-hand side tends to zero as n — > oo because &v, 
and k m is continuous. 



Theorem 2.3 Let 



m N 



3=1 



k m uniformly 
□ 



(2.62) 



where qj = ^(j — \{n N + 1)) and < iV/2. Assume that n N /N — > p 
as JV — > oo. T/ien m N ^~dq on [— 7rp, 7rp] and ra n — > m, where m is 
absolutely continuous with respect to the Lebesgue measure and symmetric, 
and there exists Q G [0, tt — fj] such that supp(m) = [-Q, Q\. 
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Proof. Let F E C([-n/2, vr/2]). Then 

/n?)^)4E f (^-^ +1 )))- H P F(q) ^ (2 - 63) 

It follows that friN — > m, and we must show that m is absolutely continu- 
ous and even. The latter follows from the fact that k m is even, which is a 
consequence of the uniqueness. To prove the absolute continuity, let e > 0. 
We must show that there exists 5 > such that m(k — 5, k + 5) < e for all 
k . Now, m(k + 5,k + 5) = j k -i {ko ^ ko+s) and we have seen that k m is 
continuous and increasing: [— 7r/2,7r/2] — > [—n + /i, n — fi]. Therefore /c" 1 is 
continuous and Ve > 35 > : k^{ko — S,k + 5) C (go — ne, <?o + ire) where 
k m {ko) = <?o- Hence J k -x {ko _ SM+5) g < e. □ 
Writing the free energy (jl.27p in the form 



mp) - y 



lim minjex-- / hiL(e lk{q) ) m N (dq), 

VT-»00 I fj J 

In M{e lk{q) )m N {dq) 



1 

£2 



P 



lim min < e x — - / \nL(e lk ) m N (dk), 

I P J 

e 2 -i J\nM(e lk )rh N (dk)\ (2.64) 



we obtain 

f(P,p) = minjex-^g / lnL(e^) dg, e 2 - ^ | lnM(e^)dgj 
= min|ei-i J In L(e tk ) m(dk) , e 2 - ^ ^ lnM(e ifc ) m(dfc)J. 



(2.65) 



By transformation to the variable y = g(q) this becomes 

{1 f yi 



-2/1 

1 



C2 - 



2tt/? 



In M(e iK M)R(y)dy \. (2.66) 

2/1 
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Here y\ = g{j/2) and R{y) = g'(q) 1 is given by (I2.58[) : 

-7r/2 



r /2 da 
R(y) = K\y) + / J{y - g{q))^ 

si n I 9 / / ^ 

_R(a)<ia 



t/2 

sin(/x) 1 f yi sin(2//) 



cosh(y) — cos(/i) 27r cosh(y — at) — cos(2/i) 

(2.67) 

If we assume that y\ = +oo then this can be evaluated by Fourier transfor- 
mation as in IK] . With 



i r°° 

R( x ) = — / R( a )e iax da (2.68) 
2vr J_ 00 



we have -R(tr) = e ^ — e 2tlx R(x) and hence 

Mx) = — r. (2.69) 

v ' 2cosh( / ux) v ; 

This is consistent with the fact that the minimum in (12.651) is attained at 
p = \ since / is convex in p. Indeed, 



1 f Vl 

R(0) = — / R(a)da = m = p. 
2vr J_ yi 



Now, R(y) = (g r Y(y) determines g 1 given that g 1 (0) = 0. By uniqueness 
of g(q), it must be the solution. Eventually, one finds 

Corollary 2.2 Assume A G [0,1). Then, the free energy is given by 

^ 1 f°° sinh[(// + w)x] sinh[(-7T — p)x] ^ 
(3 J _ OQ 2xsinh[7rx] cosh[/ix] 
= e 2 - \ ^h[ (/ ,-^] sinh[(vr-^] ^ 
p J_ 00 2xsmh[nx\ cosh[/ixJ 

where the parameter w is defined by 

a : b : c = sin|(/i — w) : sin |(yU + w) : sin(yu). (2-71) 
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3 The case A <C -1. 

For A < — 1 we write A = — cosh A, assuming A > 0. Define the new function 
a(q) by 

A _ „—ia(q) 

e ik {g) = e 

e \-ia(q) _ I ^ ' 

Then k(q) = u)\(a(q)) where u\ : [— ir, ir] — > [— n, vr] is an increasing function 
given by 



LOxiX 



f x sinh(A) , / tan (x 2) \ , . 

/ — i 7T^ / ^ du = 2 tan" 1 ,\ / /. . 3.2 

J cosh A - cos (u) \tanh(A/2) / 



In terms of a(q) the Bethe Ansatz equations read 

LO\(a(q)) = q+ u 2 \{a(q) - a(q')) m(dq') . (3.3) 

J-tt/2 

Here the function uj 2 \ is defined as in (13. 2p with the understanding that for 
| x | > it the integral expression is assumed so that u 2 \ is continuous. The 
measure m G Ai b + [— f , |] satisfies | |m| | < |. For the finite lattice, it is given 
by 

1 » n+l-2j 
m= A7Z^^' <?i = ^ 7T- (3-4) 

3=1 

Theorem 3.1 Assume A > Ao, where Ao = ln(3 + 2y / 5), w/ien ||m|| < 1/2. 
Then, for any measure m G Ai b + [— |, f] ||m|| < | ; £/iere exists a 

unique function a G L°°(m) suc/i i/iai ( fff.ffi) /ioWs /or a// q G supp(m). 
Moreover, if m is symmetric then a extends uniquely to a continuous func- 
tion on [— 7r/2,7r/2] with values in [— ir, ir] which satisfies $3.3\) for all q G 
[-7r/2,7r/2]. 

Proof. We expand tan -1 (tan(x/2)/ tanh(A/2)) into a Fourier series. For 
this, we compute first 



sinh(A) 



2ir J_ w cosh(A) — cos(x 

Hence 



- e inx dx = e" A|n| 



/ S ! nh(A) . . = Y e-We-*™ = 1 + 2 y e~ Xn cos(nx) 
cosh(A) — cos(x) ^— J 
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Integrating, we obtain the following Fourier expansion: 

00 e -An 

u\(x)=x + 2) sin(nx). (3.5) 

' n 

n=l 

Inserting this into the Bethe Ansatz equation ( 13. 3ft . we have 

f \ q 1 f / t\ / t /\ 2 ^ „ A sin(na(g)) 

a « = 1 11 11 " 1 O / oc{q')m{dq') - - ^ V e"" A + 

1 — ||m|| 1 — ||m||y I - ll m ll ~f n 



+ Ve- M " / sin(n( a (^-«( g Q)) mW) 



n=l 



(3.6) 



Notice next that it follows from (13.61): 



J a(q')m(dq') = j q'm(dq') - 2^T^ J e~ nX sm(na{q'))m (dq'). 



Let us introduce an operator: 
T[a} = 



/2 00 ^ \ f 
q'midq) H n — n > - / e~ nA sm(na(q'))m(dq') + — 
1 — I \m\\ ^— ' n J 1 



m / 1 — m ■< — ' n / 1 — m 

M M n= i ^ mi 

2 y, c _ wA sin(na(g)) | 2 y^ c -2A» f sm(n(a(q) - a(g'))) ^ 

1 — llmll ' n 1 — llmll ' / n 

11 11 n=1 11 11 n=1 </ 

We consider ( 13.61) as a fixed point problem a = T[a] and show that the map 
a 1— > T"[a] is contraction w.r.t. the L°° norm for sufficiently large A. This is 
straightforward: 



1 + 1 


m 


1 2 


1 - 1 


m 


| e A - 1 



4|| 


m\\ 


1 


1 - 


\m 


1 1 e 2A - 1 



||T(ai) -T(a 2 )|| < z rn^3 7 + ^ rn^33; 7 IK -"2 



Clearly, the factor in front of | \a\ — 0:2 1 1 is less than 1 if A > ln(3 + 2\/5). 

The same inequalitie holds for functions a G C([— f , |]), so there also ex- 
ists a unique continuous function satisfying ( 13. 3ft for all g G [— |, |]. Clearly, 
restricting this function to supp(m) yields the solution a G L°°(m). Finally, 
note that if m is symmetric, uniqueness implies that the function a must be 
odd. This in turn implies that a(q) G [— 7T, 7r] for if ai 02 G [0, 7r] then 

^2a("1 - «2) + ^2a(«1 + "2) < 2^2a("i) 

as follows easily by differentiation. □ 
We now also have an analogue of Theorem 12.31 
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Lemma 3.1 The map m h- > a m defined by Theorem \3.1\ is continuous. 
Proof. This follows from (13.61) using 

sin(na mi (q'))mi(dq') - / sm(na m2 (q'))m 2 (dq')\, 



cos(na mi (q'))m 1 (dq') - J cos(na m2 (q'))m 2 (dq')\ < ^n\\mi - m 2 \\. 

□ 

Corollary 3.1 If m n — > m weakly, and rh n is the image measure of m n 
under the mapping k mn then m n — > m = k m (m). 

As in the case A £ [0, 1) we can now conclude that the free energy is 
given by 

/CM = minjex-^g J "in L{e ik ^)dq, 

C2 "2^/ lnM ( e *)^}" ( 3 - 7 ) 
Transforming to the variable a we have 

f(e 1 ,e 2 ,e 3 -(3) = mm ^ - ^ J\nL(e^)R(a)da, 

e 2 --i-/ mM(e iu> ^)R{a)da\, (3.8) 

where R(a) = a'(g) -1 . Again, it can be evaluated by Fourier transformation, 
but now on [— n, 7r]: 

R P = T7^- (3-9) 

1 2cosh(Ap) v ; 

The resulting free energy is 

Corollary 3.2 Assume A = — cosh(A) with A > Ao- Then, the free energy 
of the 6-vertex model with periodic boundary conditions is given by 

( 1 \ \ + v ^ smh.[(\ + v)p]e- p X 

}{P > ~ ei ~P\~lT + pcosh[p\] 

~ ^--p\— + 2- i P cosh[pA] /' (3 ' 10) 

where the parameter v is given by 

a : b : c = sinh|(A — v) : sinh±(A + t>) : sinh(A). (3-11) 
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4 Concluding remarks. 



For values of A < —1 which are not large negative, numerical iteration of 
the equations f l 3 . 3 H with m = itln seems to indicate that there is in fact a 
unique solution. We have so far not been able to prove this, although it is 
possible to show that the solution is unique and increasing for small \q\. 
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